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ABSTRACT 

We have measured half-light radii, Vh, for thousands of globular clusters (GCs) belonging to the 
one hundred early- type galaxies observed in the ACS Virgo Cluster Survey and the elliptical galaxy 
NGC 4697. An analysis of the dependencies of the measured half-light radii on both the properties 
of the GCs themselves and their host galaxies reveals that, in analogy with GCs in the Galaxy but 
in a milder fashion, the average half-light radius increases with increasing galactocentric distance or, 
alternatively, with decreasing galaxy surface brightness. For the first time, we find that the average 
half-light radius decreases with the host galaxy color. We also show that there is no evidence for a 
variation of with the luminosity of the GCs. Finally, we find in agreement with previous observations 
that the average rh depends on the color of GCs, with red GCs being ~ 17% smaller than their blue 
counterparts. We show that this difference is probably a consequence of an intrinsic mechanism, 
rather than projection effects, and that it is in good agreement with the mechanism proposed in 
Jordan (2004). We discuss these findings in light of two simple pictures for the origin of the rh of GCs 
and show that both lead to a behavior in rough agreement with the observations. After accounting for 
the dependencies on galaxy color, galactocentric radius and underlying surface brightness, we show 
that the average GC half-light radii (rh) can be successfully used as a standard ruler for distance 
estimation. We outline the methodology, provide a calibration for its use, and discuss the prospects 
for this distance estimator with future observing facilities. We find (rh) = 2.7 ± 0.35 pc for GCs with 
{g — z) — 1.2 mag in a galaxy with color (g — z)gai = 1.5 mag and at an underlying surface z-band 
brightness of /i^ = 21 mag arcsec^^. Using this technique, we place an upper limit of 3.4 Mpc on the 
l-fj line-of-sight depth of the Virgo Cluster. Finally, we examine the form of the r/i distribution for 
our sample galaxies and provide an analytic expression which successfully describes this distribution. 



Subject headings: galaxies: elliptical and lenticular, cD 
general 



galaxies: star clusters — globular clusters: 



1. INTRODUCTION 

The half-light radius, r^, represents one of the most 
fundamental structural properties of globular clusters 
(GCs). Numerical simulations have shown this quantity 
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to be fairly constant throughout the dynamical evolution 
these stellar systems (Spitzer & Thuan 1972; Lightman 
& Shapiro 1978; Murphy, Cohn & Hut 1990; Aarseth & 
Heggie 1998). Moreover, there are reasons to believe that 
the present day half-light radii are faithful records of the 
characteristic sizes of the proto-GC clouds; thus, mea- 
surements of half-light radii provide unique constraints 
on the properties of the GCs at the time of their forma- 
tion (Murray & Lin 1992). 

Until recently, measurements of GC half-light radii 
were confined to Local Group galaxies where the GCs 
could be resolved in ground-based images. In the Milky 
Way, it has long been known that r^ increases systemat- 
ically with galactocentric distance, Rgc (van den Bergh 
1956). van den Bergh, Morbey & Pazder (1991) found 
a scaling relation of the form r^ oc R^-^ based on an 
analysis of 98 Galactic GCs. This situation changed dra- 
matically with the launch of the Hubble Space Telescope 
(HST) which can partially resolve the spatial profiles of 
GCs well beyond the Local Group; indeed, HST opened 
the study of GC sizes in external galaxies out to distances 
of ~ 30 Mpc (e.g., Kundu et al. 1999; Puzia et al. 1999; 
Larsen et al. 2001). An interesting result to emerge from 
these studies is the realization that the half-light radii of 
metal-rich GCs in luminous early-type galaxies appear to 
be ~ 20% smaller than those of their metal-poor counter- 
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parts — a finding that has alternatively been interpreted 
as the result of projection effects (Larsen & Brodie 2003) 
or the combined effects of mass-segregation and the de- 
pendence of stellar lifetimes on metallicity under the as- 
sumption that their half-mass radius is constant (Jordan 
2004). Unfortunately, the fields of view have been typ- 
ically too small to explore possible dependencies of 
on projected Rgc, only in one galaxy (NGC 4365) has a 
clear increase of r/j with projected Rgc been found in the 
outermost regions (Larsen et al. 2001). 

The prospects for using the sizes of GCs for distance 
estimation dates back at least to the work of Shapley & 
Sawyer (1927), who state that "measures of diameter [of 
GCs] are chiefly useful in the determination of relative 
distances". Recently, Kundu & Whitmore (2001) used 
archival HST/WFPC2 observations for a sample of 28 
galaxies to argue that the median half-light radii for GCs 
in early-type galaxies is constant at (r/j) = 2.36 ±0.4 pc. 
This led them to suggest that (rh) could be used as a 
distance indicator, although they noted that their sam- 
ple showed a tendency for the physical size of GCs to 
increase with galaxy distance, a trend they attributed to 
instrumental effects. 

In this paper, we use the GCs from our large and homo- 
geneous survey of early-type galaxies in the Virgo Cluster 
to explore the dependencies of half-light radii on various 
properties of the GCs and their host galaxies. Wc de- 
scribe the observations and data reductions in §2. An 
exploration of the dependencies of half-light radii on the 
properties of the GCs and their host galaxies is given in 
§3. Theoretical implications of the measured half-light 
radii are briefly discussed in §4. After accounting for 
various environmental dependencies, we show in §5 that 
(rh) can be successfully used as a distance indicator and 
we provide a method and calibration for its use. In §6, we 
study the form of the observed rh distribution in early- 
type galaxies and present an analytic expression which 
successfully describes this distribution across our sample. 
We summarize and conclude in §7. An appendix gives a 
description of the algorithm used to measure structural 
parameters, including half-light radii, for GCs in the ACS 
Virgo Cluster Survey. 

2. OBSERVATIONS 

A sample of 100 early-type galaxies in the Virgo cluster 
was observed as part of the ACS Virgo Cluster Survey 
(ACSVCS; Cote et al. 2004). Each galaxy was observed 
in the F475W (~ Sloan g) and F850LP (~ Sloan z) band- 
passes for a total of 720 s and 1210 s respectively. Ad- 
ditionally, the X-ray faint early-type galaxy NGC 4697 
was observed with an identical strategy as part of a joint 
Chandra- HS T progvam (GO-10003; PI C.L. Sarazin). 
Data reductions were performed as described in Jordan 
et al. (2004a). For all GC candidates, a half-light radius 
is measured by fitting PSF-convolved King (1966) mod- 
els to the observed light distribution. The methodology 
used to measure r^, is described in Appendix A. Profile 
fitting is carried out independently in both the F475W 
and F850LP filters. For this work, wc estimate r^ by 
averaging the results measured in the two bandpasses. 
Having two independent measurements in different bands 
lets us estimate systematic uncertainties arising from the 
adopted PSF models. We find that the measurements in 
the two bands have systematic differences of ^ 0.05 WFC 



pixels, or ~ 0.2 pc at the distance of Virgo (Mei et al. 
2005b). We have thus added in quadrature a systematic 
error of 0.05 pixels, or 0'.'0025, to the measured r^ values. 

The construction of a clean and complete GC catalog 
is an essential but difficult part of all studies of extra- 
galactic GC systems. While the full details of the GC 
selection process will be described elsewhere, along with 
the presentation of the GC catalogs (Jordan et al., in 
preparation), we briefly outline the technique here. Mak- 
ing use of the measured rh and .z-band magnitude and 
catalogs of the expected contamination taken from blank 
control fields and constructed for each galaxy individu- 
ally, we select GCs via maximum likelihood estimation 
of a mixture-model of GCs and contaminants. Due to 
the difficulty of distinguishing extended GCs from back- 
ground galaxies, sources with rn > 0'.'128 (~ 10 pc at the 
mean Virgo distance) are eliminated as GC candidates 
before estimating the mixture model. This method as- 
signs to each candidate a probability, Pgc, that it is a GC. 
In this work, we consider only those sources for which 
Pgc > 0.5. 

When analysing the rh measurements, it is crucial to 
minimize or avoid any selection and contamination ef- 
fects. In order to create a clean sample for structural 
parameters studies, we have restricted our sample to GC 
candidates to those with 

z < 22.9 mag , . 

0.6 < (.g - z) < 1.7 mag ^ ' 

Note that the selection on magnitude corresponds 
roughly to the expected turnover of the GC luminosity 
function at the distance of Virgo. Additionally, we have 
imposed on each program galaxy a selection on projected 
galactocentric distance. An inner cut, determined with 
the help of our artificial objects test, ensures that the 
underlying surface brightness of the galaxy is such that 
our completeness for objects with rh = 0'.'128 is greater 
than 90%. An outer cut, determined using our control 
fields, ensures that the expected contamination in our 
GC samples from background galaxies is at most 5% of 
the objects in the sample. The average expected contam- 
ination for the samples is ~ 3%, with the values increas- 
ing with decreasing luminosity: i.e., from ~ 1% for the 
giants to ~ 5% for the least luminous dwarfs. For our 
analysis, we considered only those galaxies whose final 
GC samples contained five or more objects.^'' This left 
a final sample of 67 galaxies including NGC 4697. The 
typical random uncertainties for the measured half-light 
radii are « 0'.'003 w 0.24 pc. 

As part of the ACS Virgo Cluster Survey, wc have mea- 
sured distances for most of our target galaxies using the 
method of surface brightness fluctuations (SBF; Tonry & 
Schneider 1988). The reduction procedures for SBF mea- 
surements, feasibility simulations for our observing con- 
figuration, and the calibration of the z-band SBF method 
have been presented in Mei et al. (2005ab). The com- 
plete catalog of SBF distances will be presented in Mei 
et al. (2005c). With distances to most our targets in 
hand, we may convert the observed half-light radii from 
angular into physical units, and thereby examine directly 

Additionally, four faint galaxies located close to luminous gi- 
ants were eliminated since their own GC systems are overwhelmed 
by those of their companions. These galaxies are VCC 1297, 
VCC 1199, VCC 1192 and VCC 1327. 
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the physical sizes of the GCs. In what foUows, we con- 
sider only those GCs which belong to the galaxies with 
measured SBF distances. Our final sample is summa- 
rized in Tabled which lists the 67 galaxies upon which 
our analysis is based, the blue magnitude, Bg, from the 
Virgo Cluster Catalog of Binggeli, Sandage & Tammann 
(1985), the number of GCs satisfying the above selection 
criteria, the mean GC half-light radius in arcseconds and 
its uncertainty (not including systematic uncertainties) 
and alternative names for our sample galaxies. 

3. ENVIRONMENTAL DEPENDENCIES 

As mentioned in §1, GC half-light radii are believed to 
reflect the physical conditions of the proto-GC clouds at 
the time of their formation. Thus, any systematic varia- 
tions in half-light radius can provide potentially powerful 
insights into the mechanism(s) of GC and galaxy forma- 
tion. 

The various factors that can, in principle, affect the 
half-light radii of GCs can be naturally divided into fac- 
tors that change rh individually in each GC, factors that 
affect the GC system within a galaxy, and factors that 
change it across galaxies. In what follows, we shall re- 
fer to these factors as internal, local and global factors. 
The variables we will use to parameterize internal effects 
are the GC colors, {g — z)gc, and z-band magnitudes. 
Local factors will be parameterized using the variables 
log(r/re), where r^ is the effective radius obtained from 
our isophotal analysis (Ferrarese et al. 2005, in prepa- 
ration), and the underlying z-band surface brightness of 
the galaxy, /Zz, at the position of each GC. Note that /Zz 
is preferred over its g-band counterpart because it is a 
better tracer of the underlying stellar mass and is less 
sensitive to dust obscuration. Finally, global parame- 
ters are taken to be the absolute blue magnitude of the 
host galaxy, Mb, and its corresponding color, {g — 2;)gai 
(Binggeli, Sandage & Tammann 1985; Cote et al. 2004; 
Ferrarese et al. 2005).^^ We use our SBF distance mod- 
uli (Mei et al. 2005bc) to calculate absolute magnitudes 
on a galaxy- by-galaxy basis. 

Many of these properties are highly correlated, so the 
extraction of the intrinsic dependence on each of them 
is not entirely straightforward. For example, not only 
are Mb and {g — z)gai correlated, but both quantities 
are themselves correlated with the mean color of the GC 
systems and the fraction of GCs which belong to the 
metal-rich component (with brighter and redder systems 
having a higher fraction of metal-rich GCs; e.g., Peng 
et al. 2005). Moreover, the dependence of the average 
half-hght radius, (rh), on {g — z)gc (see below) will affect 
the inferred dependence of (rh) on fiz through differing 
proportions of metal-rich and metal-poor GCs as a func- 
tion of radius: i.e., red and blue GCs are often found to 
have different spatial distributions, with the latter hav- 
ing a more extended spatial distribution (Kissler-Patig 
et al. 1997; Lee, Kim & Geisler 1998; Cote et al. 2001; 
Dirsch et al. 2003).^^ Likewise, the dependence of (r^) 

^'^ For NGC 4697, which is not a member of the Virgo Cluster 
and hence does not appear in the catalog of Binggeli, Sandage 
& Tammann (1985), we use the integrated blue magnitude from 
NED: B = 10.14 mag. 

While we will refer throughout this work to (r^) as the average 
or median of r/i, we actually calculate its value using the biweight 
location estimator (Beers et al. 1990). 



TABLE 1 

Mean Half-Light Radii for ACS Virgo Cluster Survey 

GALAXIES AND NGC 4697 
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1226 


9.31 


392 


0.033 ± 0.001 


N4472, M49 


1316 


9.58 


907 


0.032 ± 0.000 


N4486, M87 


1978 


9.81 


398 


0.029 ± 0.000 


N4649, M60 


881 


10.06 


199 


0.036 ± 0.001 


N4406, M86 


798 


10.09 


211 


0.035 ± 0.001 


N4382, M85 


763 


10.26 


257 


0.030 ± 0.001 


N4374, M84 


731 


10.51 


333 


0.024 ± 0.000 


N4365 


1903 


10.76 


190 


0.033 ± 0.001 


N4621, M59 


1632 


10.78 


244 


0.031 ± 0.001 


N4552, M89 


1231 


11.10 


136 


0.034 ± 0.001 


N4473 


1154 


11.37 


93 


0.031 ± 0.001 


N4459 


1062 


11.40 


98 


0.038 ± 0.001 


N4442 


2092 


11.51 


41 


0.037 ± 0.002 


N4754 


759 


11.80 


86 


0.036 ± 0.001 


N4371 


369 


11.80 


115 


0.037 ± 0.001 


N4267 


1692 


11.82 


56 


0.038 ± 0.002 


N4570 


2000 


11.94 


116 


0.037 ± 0.001 


N4660 


1664 


12.02 


79 


0.036 ± 0.001 


N4564 


944 


12.08 


44 


0.033 ± 0.001 


N4417 


1938 


12.11 


51 


0.041 ± 0.002 


N4638 


1279 


12.15 


74 


0.039 ± 0.002 


N4478 


355 


12.41 


8 


0.036 ± 0.003 


N4262 


1619 


12.50 


17 


0.037 ± 0.002 


N4550 


1883 


12.57 


32 


0.039 ± 0.002 


N4612 


1242 


12.60 


63 


0.036 ± 0.001 


N4474 


784 


12.67 


15 


0.033 ± 0.001 


N4379 


1537 


12.70 


8 


0.032 ± 0.003 


N4528 


778 


12.72 


25 


0.030 ± 0.002 


N4377 


828 


12.84 


16 


0.032 ± 0.002 


N4387 


1321 


12.84 


14 


0.034 ± 0.002 


N4489 


1630 


12.91 


16 


0.028 ± 0.002 


N4551 


1250 


12.91 


8 


0.041 ± 0.007 


N4476 


1146 


12.93 


36 


0.038 ± 0.002 


N4458 


1303 


13.10 


21 


0.033 ± 0.001 


N4483 


1913 


13.22 


31 


0.043 ± 0.002 


N4623 


1475 


13.36 


38 


0.038 ± 0.003 


N4515 


1178 


13.37 


51 


0.040 ± 0.002 


N4464 


1283 


13.45 


22 


0.032 ± 0.002 


N4479 


1261 


13.56 


17 


0.034 ± 0.002 


N4482 


698 


13.60 


56 


0.039 ± 0.002 


N4352 


1422 


13.64 


21 


0.044 ± 0.003 


13468 


9 


13.93 


9 


0.043 ± 0.003 


13019 


1910 


14.17 


17 


0.039 ± 0.003 


1809 


856 


14.25 


21 


0.046 ± 0.003 


13328 


1 /in 
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1087 


14.31 


29 


0.042 ± 0.002 


13381 


1861 


14.37 


25 


0.045 ± 0.002 


13652 


543 


14.39 


10 


0.043 ± 0.003 


U7436 


1431 


14.51 


28 


0.036 ± 0.001 


13470 


1528 


14.51 


24 


0.039 ± 0.002 


13501 


437 


14.54 


16 


0.054 ± 0.006 


U7399A 


2019 


14.55 


10 


0.047 ± 0.004 


13735 


1545 


14.96 


19 


0.043 ± 0.003 


13509 


230 


15.20 


11 


0.039 ± 0.003 


13101 


1828 


15.33 


7 


0.060 ± 0.005 


13635 


1407 


15.49 


11 


0.044 ± 0.003 


13461 


1539 


15.68 


17 


0.043 ± 0.002 




NC;C4697 


10.14 


229 


0.051 ±0.001 





on (g — z)gc will also influence the dependence of (rh) on 
Mb, as more luminous galaxies have redder GC systems. 

Our approach will therefore be to flrst isolate the var- 
ious dependencies (i.e., internal, local and global) on en- 
vironment, and then to "correct" the observed r/j so that 
subsequent correlations are unaffected the preceeding de- 
pendencies. We will indicate the successively corrected 
half-light radii by adding primes (i.e., r'f^,r'j^) and will 
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denote the final corrected rh (in which all environmen- 
tal dependencies have been accounted for) with a hat: 
fh- Our goals in this portion of the analysis are twofold. 
First, the individual correlations hold potentially impor- 
tant clues to the processes of GC formation and evolution 
(§0)). And second, we will show that — once the various 
dependencies on GC color, galaxy color and underlying 
surface brightness are taken into account — (fh) is re- 
markably constant across our sample, and thus appears 
to be a promising "standard ruler" for distance estima- 
tion (§0). 



ft 



3.1. Internal Factors 

One of the most remarkable properties of GCs is that 
their size appear to be independent of their mass: i.e., 
studies of the Galactic GC system show that r^ oc M^^ 
(McLaughlin 2000). In fact, this relation is equivalent to 
one of the defining equations of the "fundamental plane" 
of GCs, namely, that GC binding energy scales with lu- 
minosity as Egc oc Ll^^ (McLaughhn 2000). In FigureH 
we show r/i versus z-band magnitude for our sample of 
GCs. A regresssion analysis shows that the linear coeffi- 
cient in the function r^ = bzZ + az is not significant for 
this sample, meaning that, as expected, GC radii in the 
ACSVCS are independent of luminosity. We have also 
performed fits of this same functional form for all galax- 
ies individually. The resulting average linear terms as a 
function of absolute magnitude are shown in Figure El 
which gives the weighted average slope in each magni- 
tude bin. It is apparent that, across the factor of 350 in 
luminosity spanned by our sample galaxies, there is no 
evidence for a variation in the GC size with luminosity. 
It should be noted that this statement applies only to 
those clusters brighter than the turnover of the GC lu- 
minosity function (z < 22.9) due to the construction of 
the sample. To the extent that the z-band luminosity 
reflects GC mass, this finding provides strong evidence 
that rh oc M^^. 

On the other hand, previous HST studies of bright el- 
lipticals have shown that, in the mean, the half-light radii 
of GCs, {rh), depends on their color (g—z)gc, with the red 
GCs being « 20% smaller than their blue counterparts 
(Kundu & Whitmore 1998; Kundu et al. 1999; Puzia 
et al. 1999; Larsen, Forbes & Brodie 2001; Barmby, Hol- 
land & Huchra 2002; Larsen et al. 2001). As mentioned 
above, this color dependence will have implications for 
other rh correlations: in particular, it will affect corre- 
lations between (r/,.) and AIb, {g — 2)gai, log(r/re) and 
[iz- To avoid these complications, we choose to tem- 
porarily restrict the GC sample to blue clusters only, 
with (jj — z)gc < 1.05. Transforming from color to metal- 
licity as in Jordan et al. (2004b), this cut corresponds 
to a metallicity selection of [Fc/H] <^ —1.25. Once we 
have characterized the remaining dependencies, we will 
return to address the issue of the color dependence of 
GC half-light radii. While requiring the GCs to lie in 
a restricted range of color will not, in principle, entirely 
eliminate any dependency, it will certainly reduce it to a 
very small level. We note that the models presented in 
Jordan (2004) for the variation of rh with (g — z)gc — 
which interpret the observed difference in size between 
red and blue GCs as a joint consequence of mass segre- 
gation and depedence of stellar lifetimes on metallicity — 



20 21 22 

z (AB mag) 

Fig. 1. — Half-light radius, r^, plotted as a function of 2-band 
magnitude for the combined GCs of our sample galaxies. The data 
arc consistent with (r^) being independent of magnitude. The 
dashed line marks the median value of (r^) = 2.66 pc. 



-22 -20 -IB -16 
Mg (mag) 

Fig. 2. — Average coefficient hz calculated in one magnitude bins 
for galaxies in the ACS Virgo Cluster Survey, plotted against their 
absolute blue magnitude, Mg. The hz have been obtained from a 
linear fit of the relation log = hzZ + a io the GC data for each 
galaxy. 



predict very little variation for GCs with [Fe/H] < — 1. 

3.2. Local Factors 

Using our restricted sample of blue GCs, we now search 
for correlations between the observed rh with projected 
galactocentric radius, r^, and local galaxy surface bright- 
ness, fiz- In the two left columns of Figure |31 we show 
logr^i versus rp in units of the effective radius r^. for 
the eight brightest galaxies in our sample; in the two 
right columns, we show logr^ versus fj,z for these same 
eight galaxies. It is clear that the increases with 
both quantities, albeit rather mildly. We have performed 
linear fits of the form logr^ = a,- -I- brlog{r/re) and 
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Fig. 3. — {Two Leftmost Columns) Logarithm of GC half-light 
radius plotted as a function of log(r'p/r'e), where re is the galaxy 
effective radius measured from the ACS/WFC images. Only GCs 
with [g — z)gc < 1.05 and belonging to the eight brightest galaxies 
are shown in this figure. The dashed lines indicate the best fit- 
linear relations. Galaxy names are indicated in each panel. {Two 
Rightmost Columns) Logarithm of GC half-light radius plotted as 
a function of the galaxy surface brightness in the z-band, /x^ . Only 
GCs with {g — 2)gc < 1.05 and beloginging to the eight brightest 
galaxies are shown in this figure. The dashed lines indicate the 
best fit-linear relations. Galaxy names are indicated in each panel. 



logr/i — + bfj^fiz for all galaxies in our sample. We 
show the resulting values 6^ and 6^ , plotted as a function 
of galaxy absolute magnitude, Mb, in Figure^ We con- 
clude that most galaxies in the sample show evidence for 
an increase in rh with increasing rp/ve and decreasing 
surface brightness. A galaxy luminosity term included 
in the regression analysis is not significant at 99% level 
in either case, so we conclude that a single slope suffi- 
ciently describes these increases across the entire sample. 
The weighted average slopes for the fitted relations are 
= 0.016 ± 0.003 and br = 0.07 ± 0.01, which we wiU 
take as the coefficients of a universal relation for early- 
type galaxies. 

A slope of br « 0.07 is much lower than that found for 
GCs in our Galaxy, r/j oc Rg^j where Rgc is the Galac- 
tocentric distance. However, these relations are not di- 
rectly comparable because we observe only the projected 
galactocentric radii, r^, for the GCs in our target galax- 
ies. To compare with the Galactic relation, we have per- 
formed 100 random projections of the Galactic GC sys- 
tem using the positional data from Harris (1996), sim- 
ulating what would have been observed with HST/ ACS 
if the Milky Way was located at the distance of Virgo, 
and the GC sample was selected in the same fashion as 
for the ACSVCS targets. Only GCs with My < -7.4 
mag (i.e. only the brighter half) were used in order to 
reproduce the magnitude cut of our sample. To each 

If all GCs, rather than just the blue ones, are used in the 
analysis, the coefficients change only slightly, to fep = 0.017 and 
br = 0.09. The direction of the change is as expected from the 
combination of smaller (r^) and less extended spatial distribution 
for redder GCs. 



-33 -30 -18 -16 I -33 -30 -16 -16 

Mg (mag) Mg (mag) 

Fig. 4. — {Left) Average coefficient 6^ calculated in one mag- 
nitude bins for galaxies in the ACS Virgo Cluster Survey, plotted 
against their absolute blue magnitude, Mb- The bz have been 
obtained from a linear fit of the relation logr^ = bf^fi^ + a to 
the GC system of each galaxy. {Right) Average coefficient br 
calculated in one magnitude bins for galaxies in the ACS Virgo 
Cluster Survey, plotted against their absolute blue magnitude, 
Mb- The bz have been obtained from a linear fit of the relation 
log r ft = br log{rp/re) -f a to the GC system of each galaxy. 



simulated dataset we have fit a linear relation as above 
to find br, and recorded the result. These simulations 
shows that our Galaxy would have been measured to 
have br_MW ~ 0-3, with a dispersion around this value of 
~ 0.15. Thus, the slope we find for GCs in early-types 
in Virgo is significantly shallower than that for GCs in 
the Galaxy. While we will not explore this issue further 
in this work, it is an interesting discrepancy that may 
point to a difference in the formation or dynamical evo- 
lution of GC systems in early- versus late-type galaxies, 
or perhaps, in group versus cluster environments. 

While either log (r^ /re) and may be use to pa- 
rameterize the outward decline in rh, it is clear that 
they are not independent variables: i.e., quite generally, 
fi oc logF{rp/re), where F is an arbitrary continuous 
function with finite total integral. This does not make 
the descriptions in terms of the two variables necessar- 
ily equivalent but it renders them correlated and forces 
the choice of one. We choose to work with as the in- 
dependent variable since it avoids additional uncertain- 
ties that arise due to errors in the determination of r^, 
and in the assumption of spherical symmetry (which is 
a poor assumption for a number of our program galax- 
ies). Moreover, the use of /i^ accounts implicitly for the 
non- homologous nature of the galaxy profiles (e.g., Caon, 
Capaccioli & D'Onofrio 1993), a behavior which a simple 
power-law cannot capture. 

Using the relation between rh and fj.z, we correct the 
value of rh measured for each GC to the value ex- 
pected at an underlying galaxy surface brightness of 
Hz — 21 mag arcsec"^. The resulting corrected half-light 
radius is given then by 

r;, = rU0-°-°i6(^--2i). (2) 

Note that the form of the correction tacitly assumes 
that only the mean of the distribution is changing 
according to the derived relation. 

3.3. Global Factors 

Having examined the internal variation in rh, we now 
turn our attention to global dependencies. In this case, 
our sample is again restricted in GC color, and the half- 
light radii have been corrected as given in Equation |5J 
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(g-z),„| Mj (mag) 

Fig. 5. — (Left) Mean GC half-light radius, (r^), corrected to 
an underlying surface brightness of fi^ = 21 mag arcsec"^, plotted 
against galaxy color, (g — 2:)gai. The dashed line is the least-squares 
line of best fit, showing a residual dependence of (rj^) on (g — ^)gai- 
(Right) (rj^) as a function of galaxy absolute blue magnitude, Mg. 
The dashed line is the least-squares linear fit to the data. The dot- 
dashed curve is the relation obtained by folding the linear relation 
seen in the left panel with a quadratic relation between (g — 2)gai 
and Mg inferred from our sample. 



i.e., we are working with rj^ values. Any remaining de- 
pendency therefore must arise from global factors. In 
FigureEl we show the average half light radius (rj^) plot- 
ted against galaxy color {g — 2;)gai and galaxy blue lumi- 
nosity Mb (left and right panels, respectively). In both 
cases, there may be residual trends, particularly among 
the bluer and fainter galaxies, which seem to have some- 
what larger (rj^). Linear fits of log(rjj) with galaxy color 
and absolute magnitude give 

log(r;^) = (0.459±0.008)-(0.167±0.054)[(5-z)gai-1.5] 

(3) 

and 

log(r;j ^ (0.465±0.010) + (0.004±0.005)(Mb+20). (4) 

The coefficient of Mb is not statistically significant, 
while the color coefficient is significant at a 3cr level. Due 
to the well-known correlation between galaxy color and 
luminosity (Baum 1959; van den Bergh 1975), the use of 
either Mb or {g — z)gia should be roughly equivalent as 
a tracer of the global variation in rj^. For the ACSVCS 
sample, we can determine the relationship between these 
parameters explicitly: FigureElshows (g — 2)gai as a func- 
tion of Mb, along with the best-fit quadratic relation 

{g - z)gai = -3.2 - 0.45Mb - O.OlOAfI, (5) 

which is valid in the range —22 <^ Mb <^ —15. The need 
for a quadratic term is a reflection of the fact that the 
relation flattens out at the high galaxy luminosities, with 
a plateau in color after a certain luminosity threshold is 
crossed (cf. Tremonti et al. 2004). This might provide 
a clue as to why log(rjj) shows little or no dependence 
on an Mb term as compared with (g — z)gai. In any 
case, due to the better behavior and the added benefit of 
being independent of distance, we will choose to model 
this mild dependence using the galaxy color as opposed 
to its blue luminosity. 

We now use the dependence on galaxy color to further 
correct the half-light radii: 

r',[ = j,^10-"-016(^--21)+0.17[(s-z)gai-1.5]_ (-g^ 



CO 




-22 -20 -18 -16 
Mg (mag) 

Fig. 6. — Galaxy color, (g— z)gaii versus absolute blue magnitude, 
Mb I for galaxies in our sample. The dashed line shows the best-fit 
quadratic relation in the range —22 <^ Mb ^ —15. 



3.4. Dependence on GC Color 

Using the half-light radii corrected according to Equa- 
tion El we now return to the question of the color de- 
pendence of half-light radii. In doing so, we drop the 
restriction to blue GCs imposed in § 13.11 

To quantify variations of r^' with color, we must have 
enough GCs in any given color bin to get an accurate 
measure of the mean behavior. This requirement imposes 
a limitation on the luminosity of the galaxies to be stud- 
ied, as less luminous galaxies have GC systems which are 
almost exclusively metal-poor (i.e., many of the dwarfs 
have only a handful of GCs with {g — z)gc ^ 1.2). Thus, 
we have decided to create composite GC sets by combin- 
ing galaxies in color bins created by demanding at least 
100 GCs with colors in the range 1.3 < (5 — z)gc < 1-5. 
In addition, we consider only those galaxies which have 
a number of GCs in the range 1.3 < (5 — z)gc < 1.6 
that is at least 10% of the number found in the range 
0.8 < {g — z)gc < l-l- The motivation for this condi- 
tion is that adding galaxies which have purely metal- 
poor GC systems adds no information to the charac- 
terization of the r'^ color dependence; indeed, because 
of uncertainties in the measured distances, such galax- 
ies would add scatter only in the blue population while 
contributing few or no red GCs. Figure [7| shows the 
resulting GC samples, which in some cases consist of a 
single galaxy. Also shown are the best-fit linear relations, 
logrjj' = a{g — z)gc + b (transformed to rj^'), the median 
values of r'/^ in bins of {g — z) along with their 99% con- 
fidence intervals, and a model from Jordan (2004) with 
central potential^^ Wq = 9, normalized such that the 
model at [Fe/H]= —1.5 coincides with the median r^' at 

Jordan (2004) modelled the effects of mass segregation using 
multimass isotropic Michie-King models (Gunn & Griffin 1979), 
which are characterized by their central potential Wq, which is the 
gravitational potential at the center of the cluster normalized by a 
characteristic velocity variance. A value of Wo = 9 is comparable 
to a value of c ~ 1.6 for the concentration of a single-mass isotropic 
King (1966) model, with the concentration increasing with Wq. 
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Fig. 7. — GC half-light radii, r^', corrected for the effects of 
surface brightness and galaxy color, plotted against GC color for 
GC subsamples formed by combining galaxies in color bins as de- 
scribed in the text. The color bins, or the name of the individual 
galaxies if only one present, are recorded in the different panels. 
The number of galaxies used to create the GC subsamples is in- 
dicated in brackets. The large symbols indicate the median rj^' in 
different GC color bins. The error bars indicate the 99% confidence 
intervals. The dashed line represents the best linear fit of the form 
logr^' oc [g — z)gc to the data (transformed to r'^). The solid line 
represents the prediction of the behavior of rj^' with (g — z)gc from 
a model from Jordan (2004) with a central potential Wo = 9. For 
the bin with {g — 2)gai < 1-52 we have included a second model 
with Wo = 7 (the upper one). 



(ff - 2)gc = 1-05 mag. 

It is apparent that the rj^' show a significant correlation 
with {g — z)gc for all of the subsamples in Figured Ex- 
cept for the group of bluer galaxies shown in the bottom 
left panel, the behavior is very consistent, with 

logr;'cx(-0.17±0.02)(g-z)ge, (7) 

valid in the range 0.8 <i {g ~ z)gc 1-6. The bluer bin 
shows a somewhat shallower behavior, with a slope of 
—0.09 ± 0.04, but still within ~ 2a of the mean behavior 
of the redder galaxies. Because the bluer bin required 
25 galaxies for its construction, it is much more prone 
to uncertainties arising from errors in the measured dis- 
tances. Futhermore, these galaxies contain few red GCs, 
so their usefulness in probing the relation between r^' and 
{g — z)gc is limited. We will therefore assume that the 
behavior exhibited by galaxies with {g — 2;)gai ^ 1.52 is 
valid universally, and will use this relation to perform the 
final correction of the GC half-light radii (see below). 

It is remarkable how well a single model from Jordan 
(2004), with Wq — 9, does in matching the data in Fig- 
ure 13 With the possible exception of the bin contain- 
ing the bluest galaxies, it reproduces the observed trends 
very well. As mentioned in only two mechanisms have 
been proposed to explain the observed tendency for red 



GCs to be smaller than their blue counterparts. Larsen 
&: Brodie (2003) proposed that this effect could be the 
consequence of projection effects, combined with the dif- 
fering spatial distributions of the red and blue GCs. The 
red GCs are observed to be more centrally concentrated 
so that, if the GCs follow a relationship between galac- 
tocentric distance and {r^) which is similar to that of 
GCs in the Milky Way, projection effects can lead to the 
red GCs appearing systematically smaller (due to the 
fact that they are, on average, physically closer to the 
galaxy center). As noted by Larsen & Brodie (2003), 
however, this mechanism requires some fine-tuning to 
work. Jordan (2004), on the other hand, has proposed 
that the observed size difference can be explained as 
the consequence of mass segregation within individual 
GCs, combined with the dependence of stellar lifetimes 
on metallicity. Simply stated, if the half-mass radii are 
assumed to be independent of metallicity, the half-light 
radii will, at fixed age, differ because the brightest stars, 
which dominate the overall light profile, will be more 
massive in red (metal-rich) GCs and thus more centrally 
concentrated. 

We can use our new measurements to discriminate be- 
tween these alternatives. First, we note that our observa- 
tions point to a significantly shallower relation between 
r/i and galactocentric radius in early-type galaxies than 
the one observed in our Galaxy. This latter relation was 
assumed by Larsen & Brodie (2003) and is necessary for 
the projection mechanism to explain the observed trend. 
Using shallower relations they found that projection ef- 
fects are unable to reproduce the observed size difference. 
Thus, our measurements seem to rule out projection ef- 
fects as the primary mechanism responsible for the ob- 
served difference, although it is certainly possible that 
they play some part in its origin. 

In addition, the importance of projection effects is ex- 
pected to depend on the fraction of the galaxy surveyed, 
as they are expected to decrease strongly beyond one 
effective radius of the GC system, and to disappear en- 
tirely as rp ^ oo (Larsen & Brodie 2003). The galaxies 
with (5 — 2;)gai '> 1.52 span an enormous range in riii„/re, 
where rum is the upper radii for inclusion of GCs in each 
galaxy and re is the effective radius of the galaxy light. 
These galaxies have 1 <^ ruin/re ^ 20, with a median of 
Trnn/fe = 3.3 and an interquartile range of ^ 8. That the 
observed difference in is observed to be nearly identi- 
cal among these galaxies argues strongly for a mechanism 
that is independent of projection effects. It might be ar- 
gued that the fact that the decrease in size seems to be 
milder for the bluest galaxies is consistent with at least 
some variation due to projection effects. But the distri- 
bution of riim/re for the galaxies in the bluest bin is not 
extremely different from those in the redder ones, hav- 
ing 2 ^ riim/re <^ 25, a median of rum/re = 4.8 and an 
interquartile range of 5.5. Thus, it would be fortuitous 
that the change of regime occurs between samples with 
median nim/^'e of ~ 3.3 and ~ 4.8. 

The mechanism of Jordan (2004) seems to fare better, 
despite the fact that it is a clear oversimplification to 
compare the data with a single model with Wq — 9. In 
reality, GCs will span a range in Wq , and there could well 
exist a systematic difference in the average central poten- 
tial of GCs belonging to the blue (low-luminosity) galax- 
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Fig. 8. — (Upper Panels) Fully corrected mean GC half-light 
radius, (r^), in parsecs, plotted as a function of galaxy color, [g — 
2)gab in thfi left panel, and against and absolute blue magnitude, 
Mfl, in the right panel. The dashed line indicates the biweight 
estimate of the mean, {fh) = 2.72 pc. The dotted lines indicate 
the dispersion of 0.35 pc. {Lower Panels) Fully corrected mean 
GC half-light radius, (r^)) ii^ arcseconds, plotted as a function of 
galaxy color, {g — ^)gai, in the left panel, and against and absolute 
blue magnitude, Mg , in the right panel. The dashed line indicates 
the biweight estimate of the mean, (f^) = 0''033. The dotted lines 
indicate the dispersion of 0"0035. VCC731, a member of the W 
Cloud, has been excluded from the computation of these quantities. 



ies. As we show in Figure 13 a model with Wq = 7 for 
the bluer bin is able to reproduce the observations quite 
well. In any case, the overall conclusion to be drawn from 
the data in Figure [3 is that, whatever is responsible for 
the size difference between red and blue GCs does, this 
difference does not seem to arise mainly from projection 
effects but is instead the result of some property intrinsic 
to the GCs. 

3.5. Combined results 

Having systematically explored the variations in half- 
light radii with a number of internal, local, and global 
factors, we now use the observed dependencies to correct 
the individual measurements to the expected value for 
a GC with color (g — z)gc = 1-2 mag, at an underlying 
surface brightness of — 21 mag arecsec"^, in a galaxy 
with (g — z)gai = 1-5 mag. In other words, we assume 
that only the mean of the distribution is changing to 
the above scaling relations, so that the corrected radii, 
f/i, are given by 

fh = 7-^10-0-016(M.-21)+0.17[(g-2:)g,i-1.5]+0.17[(g-z)-1.2] 

(8) 

We plot {ffi) for our sample galaxies, in both arcsec- 
onds and parsecs, against {g — z)gc and AIb in Figure|S| 
It is clear that {fh) is remarkably constant across our 
sample galaxies. As such, this quantity appears to be 
a "standard ruler", suitable for use as an extragalactic 
distance indicator. We shall explore this issue in more 
detail in ^ 

We emphasize that none of the dependencies discussed 
above are likely to be due to selection effects. Compact 



clusters at large galactocentric radii (or, equivalently, at 
low surface brightness) would be readily detected and 
measured in our survey. The inner radial cut imposed 
on each galaxy ensures that we are highly complete for 
extended GCs. The decrease in (rj^) with galaxy color 
could be suspected to be due to contamination preferen- 
tially affecting the small GC systems of faint, blue galax- 
ies. However, we have accurate estimates of the expected 
contamination measured directly from our control fields 
(Peng et al. 2005, in preparation) so this is unlikely to 
be the case. 

While from an observational point of view, the local 
and global factors can be separated naturally, they could 
in principle arise from a single relation which determines 
rfi- As a simple illustration, let us assume that is 
related to the surface mass density, Sgai, through the 
relation (r^) oc S^^^ = [(A/gai/Lgai)/gai]~'', where /gai 
denotes the galaxy surface brightness and (Mgai/Lgai) 
its mass-to-light ratio. Clearly, in this case 

log(r,,) oc -2;log(Mgai/Lgai) + OAxfiz- (9) 

The mass-to-light ratio will vary systematically with 
both the galaxy luminosity and color. Taking 
(^^gai/-^gai) oc Lg;j from Cappellari ct al. (2005), we 
get 

\og{rh) oc 0.12a;Mi3 + OAxfi^- (10) 

where we have used the fact that Lgai oc iq-oaMb ^ Uj^_ 
der this simple hypothesis, a value of a; ~ 0.04 would be a 
reasonable description of the observations (cf. eqs 4 and 
8). This simple example serves to show that empirical 
trends can guide our theoretical understanding of both 
GC formation and galactic structure. 

We will now explore the implications of the observed 
global behavior for models which make quantitative pre- 
dictions regarding the half-light radii of GCs and their 
dependence on environment. 

4. THEORETICAL CONSIDERATIONS 

Our finding that the half-light radii of GCs show clear 
correlations with some properties of their host galaxies 
should not obscure the remarkable fact that over a fac- 
tor of ~ 350 in galaxy luminosity, (rh) varies by only 
^ 10% (after accounting for internal and local factors). 
Clearly, one could assume that, to first order, (r/j) is con- 
stant: (rh) ~ L^^i- Thus, our results demonstrate that 
any paradigm for the formation of GCs should predict 
an exponent rj in the relation (r/j) L^^j which is very 
close to zero. In this section, we show that two possible 
hypothesis for the origin of the observed Vh of GCs do 
lead to roughly this result. It should be stressed at the 
onset, however, that the level of theoretical understand- 
ing we have on process of GC formation is sketchy at 
best. Indeed, the observations are far ahead, in the sense 
that they provide an increasingly clear, albeit sometimes 
puzzling, set of constraints to which any theory of GC 
formation must conform. 

4.1. Half-Light Radii as the Result of 
Pressure- Confined Proto-GC Clouds 

It has been suggested that the half-light radii of GCs 
reflect the structure of self gravitating proto-GC clouds 
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which, at the time of GC formation, were in pressure 
equihbrium with the surrounding gas. The structure of 
such clouds has been calculated by McLaughlin & Pu- 
dritz (1996), who determine that their limiting radius 

satisfies the relation i?ii,„ oc M'^-^Pgas^'^ (cf. equation 
2.11 in McLaughlin & Pudritz 1996), where Pgas is the 
pressure of the ambient gas. If we posit that the aver- 
age rh observed today is linearly related to i?iim, then it 

follows that (rh) oc M°-^ P^as^^ , where we assume that 
GCs of differing masses that form at approximately the 
same radius in a galaxy come from clouds with differing 
velocity dispersions. It is immediately apparent that this 
hypothesis is in strong conflict with observations due to 
the factor Af°-^. The observational evidence presented in 
§3 argues for (rh) oc To resolve this discrepancy one 
might propose that the star formation efficiency scales as 
some power of the proto-GC cloud binding energy (e.g.. 
Ashman & Zepf 2001). While such a scaling relation 
could, in principle, erase the mass dependence, there is 
no physical basis to justify such a scaling. In any case, 
we will carry on by supposing that a mechanism exists to 
render this pressure support picture viable, and examine 
what this scenario predicts for the global behavior of (rh) 
with host galaxy luminosity. 

We begin by supposing that {rh) oc Pg'as^^- If the am- 
bient pressure is distributed like the dark matter, then 
averaging over an entire galaxy gives 



Af2 



^1/4 



(11) 



where Mvir and r^ii are the virial mass and virial radius of 
the dark halo. We assume now that early-type galaxies 
all formed at similar large redshifts. The halos should 
then have the same average density, or M^w oc r^jj.. Since 
the circular velocity at the virial radius is given by w^;^ = 
GAfvir/^vir, we have then that Vvh- oc Uvir and therefore 

we expect (r/j) oc v^-^^ . 

In order to progress further, we assume now that 
the dark matter halos are described by NEW profiles 
(Navarro, Erenk & White 1995). Then the circular ve- 
locity Vc profile is 



Ve{r) 



1 ln(l -I- cx) — cxj (1 -I- ex) 



(12) 



X ln(l-Hc)-c/(l-hc) 

where x = r/r^-„ and c is the concentration parameter. 
This velocity profile implies a peak circular velocity at 
radius cx ~ 2.163, at which point the value is 



\^ \ c 

-0.2162— r^—n r « 0.2162[2.742, 

Wvir/ ln(l + c)-c/(l + c) ^ 



,0.41 



.(13) 

where the second approximation is accurate to within 
~ 5% in the range 4 < c < 20. Eurthermore, simulations 
show (Bullock et al. 2001) that c ~ M~°-^^ and thus 
^ c^-^'v^i^ ~ ^iir^- Therefore, (r/,.) is related to the 
peak circular velocity of the halo through 



{rh) oc V. 



-1/2 



-0.54 



(14) 



Einally, Ferrarese (2002) connects Vp to the stellar ve- 
locity dispersion of the central bulge by 



Fp-CTfJ, a; = 0.89 ±0.04 



(15) 



and the Faber- Jackson relation has (Tq oc l}JJ^ so ulti- 
mately 

(r,)cxL-V^ (16) 

Note that a scaling relation of the form of Eauationll6l 
can also be obtained much more simply by assuming 

that ctq oc Wvir, and as consequence t^vir oc L^J^ and 

{rh) oc Wvir^^ °^ ^s&\^ ■ 1^ both cases, we get a very 
mild dependence of {rh) with L^^i^ albeit stronger than 
that which is observed. We note that if we had assumed 
that the average ambient pressure of the proto-GC clouds 
comes from gas distributed as the present-day starlight 
and used fundamental plane scaling relations, we would 
get a similar dependence but of opposite sign. That is to 
say, in this case we would posit 



{rh) oc P^^l'^ oc 



-1/4 



(17) 



where r,. is the effective radius of the galaxy light. 
The galactic fundamental plane scalings of r^ oc L^^^y 
and M/L oc L^^^y (Ha§egan et al. 2005 and references 
therein) then give (r/j) oc ^gai^- Overall, it seems fair to 
conclude that the pressure hypothesis predicts roughly 
{rh) oc il^jjj, in accord with the observations from §3. 

4.2. Half- Light Radii as a Consequence of Tidal 
Limitation 

An alternative hypothesis is that the half-light radii of 
GCs are proportional to their tidal radii, rt- This runs 
immediately into the problem of the non-homologous na- 
ture of GCs. In fact, the King concentration parameter 
c is correlated with the GC mass through the rough re- 
lation 10'^ ^ Af*^-*, which could be the reflection of an 
original "fundamental line" (McLaughlin 2000). The rh 
in a King model is a function of c times rt so that, via 
the dependence of c on Af, it follows that a simple pro- 
portionality between rh and rt for all GCs cannot be 
appropriate. Also, there is no clear mechanism to phys- 
ically transmit the information from the tidal regions to 
determine the half-light radius. Finally, if the rt are de- 
termined by the tidal field of the host galaxy, then we 
would have {rh) oc {rt) oc M^^^ (Binney & Tremaine 
1987). Thus, we would need to identify a mechanism 
that erases the mass dependence as in the pressure hy- 
pothesis above. 

Assuming that such a mechanism exists, what would 
this hypothesis predict for the global behavior? In this 
case {rh) should be related to the mean density of the 
galaxy's dark matter halo, or 



{rh) oc (Afvir/r^ir) 



-1/3 _ 



= constant 



(18) 



Therefore, this picture would also predict a roughly con- 
stant {rh) across galaxies. It is interesting that both 
hypotheses predict roughly the observed global behav- 
ior. We conclude that the current observations cannot 
discriminate between these two scenarios. Indeed, the 
observational constraints seem to be ahead of what our 
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limited knowledge of GC formation can predict, although 
it is clear that fundamental aspects of GC formation are 
encoded in them. 

5. HALF-LIGHT RADII AS A DISTANCE INDICATOR 

In the upper panel of Figure |2| we plot (r/j) in arc- 
seconds versus inverse distance, D, for those galaxies 
in our sample with SBF distance measurements (Mei 
et al. 2005c). If {rh) is constant in physical units, then 
we expect ((r?i)/l") oc D^^ at low redshift (the dashed 
line). Although the vast majority of our sample galaxies 
are located at very nearly the same distance (i.e., the 
mean Virgo Cluster distance oi D = 16.5 ± 0.25 Mpc; 
Mei et al. 2005c), two galaxies in our sample give us 
some additional leverage to test this idea: (1) VCC 731 
(NGC 4365) is a member of the Virgo W Cloud (de Vau- 
couleurs 1961) and has an SBF distance of D = 23.3±0.6 
Mpc (Mei et al. 2005c); and (2) NGC 4697, a non- 
member of the Virgo Cluster which lies in the foreground, 
with an SBF distance of 11.3±0.5 Mpc (Tonry et al. 
2001; Mei et al. 2005c). 

The upper panel of Figure El reveals a scahng of (r/j) 
with which has roughly the expected form, albeit 
with appreciable scatter for some of the fainter Virgo 
members. Note, however, that the half-light radii plot- 
ted in this panel have not been corrected for any depen- 
dence on GC or host galaxy properties as discussed in §3. 
In the bottom panel of this figure, we plot the average 
corrected half-light radii, {fh)^ against inverse distance. 
The dashed line, (fh) oc D~^, is drawn so that (fh) has 
a value of 0'.'033 at our adopted mean distance of Virgo 
(Mei et al. 2005b). Note that line is not a fit, although it 
accurately predicts the measured (f/j) for both VCC 731 
and NGC 4697. It is clear that (fh) has the potential to 
be a powerful standard ruler for distance estimation. But 
as the upper panel of Figure El shows, even without any 
correction the directly measured (rh) can be a useful dis- 
tance indicator when restricted to the luminous galaxies 
(a class to which both VCC 731 and NGC 4697 belong). 
And trading some accuracy for simplicity, the mild na- 
ture of the required corrections means that a directly 
measured (r/j) can give a rough distance measurement 
for any early-type galaxy. 

It should be noted that at least two other galaxies 
in the ACSVCS have celestial positions close to the W 
Cloud: VCC 571 and VCC 575. While they do not have 
enough GCs remaining after the selection process de- 
scribed in §121 to be included in our sample (they are left 
with just 3 GCs each), their measured {fh) are neverthe- 
less consistent with membership in the W Cloud. More 
specifically, VCC 571 has {fh) = 0'.'023, while VCC 575 
has {"fh) = 0'.'021. An SBF distance is available only for 
VCC 575, which confirms that it is at a similar distance 
as VCC 731. While the number of GCs used to measure 
{ffi) for VCC 575 is indeed low, it is reassuring that its 
measured value is consistent with it being at the distance 
of the W Cloud (« 23 Mpc), lending additional support 
to the validity of using (fh) as a distance indicator. The 

1^ Tonry et al. (2001) list the distance as 11.7 Mpc. To compare 
this value to our new SBF distance, the zeropoint must be adjusted 
to match that used in Mei et al. {2005b). As a consequence, the 
distance modulus in Tonry et al. (2001) should be shifted by -0.06 
mag, so that D = 11.4 ± 1.6 Mpc. This is in excellent agreement 
with our new measurement of 11.3 Mpc. 
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Fig. 9. — {Upper panel) (r^,) (no corrections) versus inverse 
distance D~^. Filled symbols are galaxies with Mb < —19 mag, 
while open symbols are those with Mg > — 19 mag. The star is 
NGC 4697 and the filled circle surrounded by a square is VCC 731. 
The dashed line represents the relation (r^) oc normalized 
to be 0'.'033 at D = 16.5 Mpc, and it is not a fit. Error bars 
are included only for NGC 4697 and VCC 731 for clarity. The 
asterisk is VCC 575. {Bottom panel){ff^) versus inverse distance 
D-l. Filled symbols are galaxies with Mg < — 19 mag, while open 
symbols are those with Mg > — 19 mag. The star is NGC 4697 and 
the filled circle surrounded by a square is VCC 731. The dashed 
line represents the relation {rh) oc normalized to be 0'.'033 at 

D = 16.5 Mpc, and it is not a fit. Error bars are included only for 
NGC 4697 and VCC 731 for clarity The asterisk is VCC 575. 



position in the sky and the measured (rh) for VCC 571 
suggest that it is indeed a member of the W Cloud. 

We now put forward the following method and calibra- 
tion for the use of GC half-light radii in distance estima- 
tion. If a set of {rh,i} are measured on a given galaxy 
they should first be individually corrected for their color, 
underlying surface brightness and the host galaxy color 
by computing fh,i according to equation |S1 The sur- 
face brightness can be measured from the same frame 
where the GCs are detected and analyzed, and all of the 
quantities needed to perform the correction are distance 
independent. Then, if the median of the fhA is {fh) in 
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arcseconds, the distance in Mpc follows from 

I? = (0.552 ± 0.058) (^-^^ Mpc (19) 

where the uncertainty reflects the uncertainty in the as- 
sumed mean Virgo distance and the PSF systematics. 
While the expression above is in principle applicable to 
any early- type galaxy within the color/luminosity range 
spanned by our sample galaxies [i.e., 1.1 < (5 — ^)gai ^ 
1.6 and —22 <^ Mb ^ —15.5], in practice this technique 
would be applied almost exclusively to luminous giants 
which have populous GC systems, and where large sam- 
ples of GCs can be easily accumulated and contamina- 
tion is less of an issue. Additionally, in equation |S1 mea- 
surements of the {g — z) colors for the GC and galaxy 
are required, which would entail observing in these same 
bandpasses or transforming the color in hand to (g — z), 
which would add some uncertainty. Needless to say, this 
also means that observations in at least two bands are 
needed to perform the correction. Thus, we note that 
for luminous galaxies (i.e., those with Mb ~ —21 mag) a 
distance may be obtained without the benefit of color in- 
formation by correcting the set of rh^i with the following 
approximation to equation |S1 

rh ^ rhlO-°-°^''^'^'-^'^+"-°\ (20) 

The distance then follows by applying equation 1191 The 
additional uncertainty incurred by this approximation if 
applied to a giant elliptical can be estimated to be of 
order a few percent. Of course Mb won't be known be- 
forehand and it would not be possible to know if the ap- 
proximation is valid. In practice, this complication can 
be circumvented by applying the method only to the ~ 10 
brightest early-type members of a given galaxy cluster. 

What accuracy can be expected from this method? 
The rms deviation from a constant value (r/j) = 0'.'0333 
among our Virgo galaxies is 0'.'0034 (where VCC 731 
has been excluded). Thus, a single measurement can 
be expected to give a distance with external accuracy 
of « 11%. It is interesting to see that using our sam- 
ple of (f/i) we can already set an upper limit to the 
depth of Virgo (once again, excluding the background 
W Cloud). Assuming that (fh) is strictly constant and 
neglecting observational uncertainties, the observed dis- 
persion translates into a distance dispersion of ao ^ 1.7 
Mpc, or a 1-a line-of-sight depth of 3.4 Mpc. While this 
finding is in good agreement with several previous mea- 
surements of the depth of Virgo around M87 (e.g.. West 
& Blakeslee 2000, Neilsen & Tsvetanov 2000, Jerjen et al. 
2004) , it effectively rules out the claim of Young & Cur- 
rie (1995) that the Virgo Cluster has a very elongated 
distribution (extending from ~ 8-20 Mpc) along the line 
of sight. Thus, this approach complements previous evi- 
dence against this claim from SBF measurements (Jerjen 
et al. 2004; Mei et al. 2005c). 

5.1. Geometrical Calibration and an Independent 
Distance to Virgo 

Up to now we have relied on our SBF distance measure- 
ments to establish the physical size of GCs in Virgo and 
to provide a calibration for the use of half-light radius as 
a distance indicator. Being reliant on SBF means that 



our distance measurement method is tied to the primary 
calibrators of the SBF method itself, and in particular 
to the calibration of the Cepheids. In this section, we 
calibrate the method using the Milky Way GC system 
only, thereby providing a simple geometric calibration of 
the method which is largely independent of any other 
extragalactic distance measurement method. 

First, we note that all of the dependencies used above 
to compute corrected radii, fhj, are independent of our 
SBF distances. Thus, we keep the form of these depen- 
dencies intact and attempt derive the zeropoint for the 
physical size using the Milky Way GC system. Of course, 
this approach faces has the difficulty of estimating the pa- 
rameters needed for the corrections for the Milky Way, 
and the complication that our Galaxy is a disk system 
as opposed to the early-type galaxies targeted in the 
ACSVCS. These caveats notwithstanding, we now esti- 
mate a rough zeropoint to provide a direct calibration. 

The quantities needed to determine (fh) for Galactic 
GCs are the {g — z) colors of the GCs and the Galaxy 
itself, and the surface brightness profile in z. Instead 
of working with the latter, we choose to use log(rp/re), 
where for the effective radius we adopt the value given 
by de Vaucouleurs &; Pence (1978), r^ ~ b kpc. We 
then use the Harris (1996) catalog of Galactic GCs to 
obtain the r^j, [Fe/H] and galactocentric distances Rgc- 
In order to estimate their [g — z) colors, which are not 
directly available, we transform from [Fe/H] to [g — z) 
as in Jordan et al. (2004b) (which assume an age of 13 
Gyr). Finally, we take the integrated {g — z)gai color of 
our Galaxy to be {g — z)shc — 1.17, the value appropri- 
ate for a type Sbc galaxy according to Fukugita et al. 
(1995). With these quantities in hand, we obtain a set 
of corrected half-light {rh,i} radii using the prescriptions 
discussed above, namely 

f, ■ = rh .lO-0-"'^^°s(rp/5)+Q.17[ig-z)-1.2)+Q.17(g-z)sbc-l-5] 

(21) 

for each of 100 simulated Galactic datasets with ran- 
domly projected galactocentric distance and with GCs 
selected according to our selection procedures (cf. S I3.2|I . 

From these simulations, we find {fh) ~ 2.54 ± 0.1 pc 
for the Galactic GC system. Note that the quoted un- 
certainty does not include systematics arising from our 
assumed properties of the Galaxy. These can be roughly 
estimated to be at least 15% by varying the assumed 
Te by 20% and changing the assumed color to that of 
a galaxy of type Sab and Scd. Using the fact that the 
observed (fh) for our sample is 0'.'033, we find that the 
mean distance of Virgo using this direct calibration is 
16±2.3 Mpc. While this calibration should be viewed 
with some caution (being based on a late-type galaxy 
for which we do not have a direct measurement of its 
global characteristics), it is comforting to see that the 
distance of Virgo is in good agreement with other meth- 
ods. Moreover, because this estimate circumvents several 
of the problematic issues involved in other methods used 
to establish Virgo's distance, such as the metallicity de- 
pendence of the Period-Luminosity relation for Cepheid 
variables and the distance modulus to the Large Magel- 
lanic Cloud, it provides a largely independent confirma- 
tion of the Cepheid/SBF distances. 

6. THE FORM OF THE Rh DISTRIBUTION 
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So far we have examined only the dependencies and 
properties of a location parameter of the distribution of 
GC half-light radii in early-type galaxies. In fact, in 
applying the various corrections to the set of rh.i, we 
have tacitly assumed that only (vh) was changing and 
that the shape of the distribution is invariant. We now 
show that this assumption is well justified and present an 
analytic form which accurately describes the observed 
distributions. Throughout this section, we will work with 
the fully corrected half-light radii, f/j and we additionally 
include galaxies that have less than 6 GCs after the cuts 
described in § El 

To study the form of the distribution of half-light radii 
across the color and luminosity range of our galaxies, we 
chose to create bins in galaxy color. To minimize sta- 
tistical noise, we require that each subsample contains 
at least 500 GCs (the reddest galaxy, VCC 1978, is as- 
signed a bin of its own although it does not satisfy this 
condition). The resulting fh distributions are shown in 
Figure [TUI along with the distribution of the combined 
sample. Each panel shows a Gaussian kernel density es- 
timate of the observed log fh distribution, with a band- 
width h = 0.035 (solid curves). The corresponding 99% 
confidence bands are indicated by the dotted curves. It 
is apparent that the distribution of log is very similar 
across the galaxy color range of our sample. The primary 
features are a prominent peak at r/j ~ 2.5 pc and an ex- 
tended, large-f/i tail which is much heavier than if the 
distribution was normal. The distribution is noticeably 
skewed, with an excess on the right side of the peak. 

We emphasize that the distributions shown in Fig- 
ure El are rather uncertain for fh ^ 1 pc, as the mea- 
surements in that regime begin to suffer from biases 
to slightly larger r^, due to the fact that the GCs are 
close to being unresolved. This difficulty notwithstand- 
ing, the decline toward small is certainly real. Un- 
fortunately, the data do not allow us to characterize the 
regime f/j ^ 1 pc accurately, so the following discussion 
should be taken to apply to the range ^ <^ f/i <^ 10 pc 
only. Note that the upper radius is imposed by our sam- 
ple selection; in practice, this limit is set by the difficulty 
of distinguishing extended GCs from background galax- 
ies. Thus, highly extended GCs, which certainly exist in 
our Galaxy, are not considered here. In any case, there 
appears to be no sudden upturn in the distribution for 
fft ^ 1 pc or fh ^ 10 pc (see, e.g.. Figure 12 of IIa§egan 
et al. 2005), so it is clear that the bulk of the distribution 
lies in the regime accesible to us. 

We have devised an analytical distribution function 
with four parameters which is able to describe the ob- 
served distribution of log f/j , such that the probability of 
finding a GC with logf?i in an interval of width d(logf;i) 
around f/j is given by J'(\ogfh)d{logfh)- The form of 
this function is 

j{\ogf^\e) = fg {log f^\^i,|3,)+{l-f)g{-\ogh]-^^, 1^2 

(9.9. 

where 

gix\^l,P) = iexp[-(a;-Ai)//3-exp(-(x-Ai)//3)], (23) 
e = ifi, /?!, /?2, /), /? > and < / < 119. 

The function Q can be derived as the distribution of extreme 



(g-z)g„i < [41] ::i.41 < (g-z)g„| < K52 [17] 

N=504 - N=563 



(22) 




VCC 1316 = M87 1.52 < (g-zj^^, < 1.54 [6] 
N=766 - N=992 




Fig. 10. — Distribution of log for GC samples built by combin- 
ing galaxies in galaxy color bins as described in the text. The color 
bins (or the name of the individual galaxies if only one is present) 
are indicated in the different panels. The number of galaxies used 
in creating the combined GC sample is indicated in brackets and 
the number of GCs included is indicated by the value of N. The 
bottom right panel is built by combining all galaxies in the sam- 
ple. The solid curves show Gaussian kernel density estimates with 
a bandwidth h = 0.035; the dotted curves show their 99% confi- 
dence bands. Th e thick dashed line represents the best-fit function 
J7 (Equation convolved with a Gaussian with cr = h. 



The equivalent form of this distribution in terms of fh, 
such that the probability of finding a GC with fh in an 
interval df^ around fh is given by J{fh) dfh, is 



J{fh\Q) 

where 



/g(f„|/i,/3i) + (l-/)f-^e(f-i|A-\/32) (24) 



xf3 



-1/13 



exp(~(a;/M)-i/^), (25) 



values for independent identically distributed random variables. 
While this fact might hold some clue as to the physical basis for 
the observed distribution, we use it here strictly as a statistical 
modelling function without an a priori physical justification. 
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a;,/i,/3 > 0, < / < 1 and the parameters 9 = 
(/i, (3i,P2, f) are related to the ones in J' by the relations 
p. = [3i^2 = ln(10)/3i,2 and f ^ f. These functions 
satisfy J'(logr;,) d(logr^) = J{rh) dru = I. 

We have obtained the best-fit parameters /x, /32 and 
/ via maximum-likelihood for the various subsamples 
shown in Figure ^| as well as for the combined sample 
of GCs. The resulting functions, convolved with a Gaus- 
sian with a = 0.035 in order to faithfully compare them 
with the density estimates, are shown in Figure [TUl with 
a thick dashed curve. The function is able to produce a 
close match to the observed distributions. 

The maximum-likelihood parameters obtained for each 
galaxy subsample are fairly consistent with the parame- 
ters inferred for the combined sample of all galaxies. This 
is true despite the fact that uncertainties arising from the 
effect of distance uncertainties in combining the val- 
ues for several galaxies are not included in the formal er- 
ror estimates. The parameters for the combined sample, 
with uncertainties reflecting the dispersion of the values 
inferred for all galaxy color bins, are 



and 



Pi 

f 



f 



0.407 ±0.013 
0.117 ±0.005 
0.078 ±0.011 
0.7 ±0.1 

2.55 ±0.08 
0.27 ±0.012 
0.18 ±0.025 
0.7 ±0.1 



(26) 



(27) 



We conclude that the observed distribution of fh is well 
represented by single distributions of the form 

J(logf;,.|^ = 0.407, /3i 0.117, /32 - 0.078,/ = 0.7) 

(28) 

or, equivalently, 

J{fh\fl = 2.55, ;3i = 0.27, = 0.18,/ = 0.7). (29) 

We now examine the applicability of the proposed func- 
tion to the case of the Galactic GC system. In Fig- 
ure ^2 we show a Gaussian kernel density estimate with 
h = 0.06 of the distribution of logf^ for all GCs in our 
Galaxy satisfying fh < 10 pc and My < —7.4 mag, along 
with the 99% confidence bands. The individual r^^i have 
been corrected as described in §0 We also show the best- 
fit form of i7 shown in equation [2H1 and convolved with a 
Gaussian with a = 0.06 as a dashed curve. Without any 
adjustment the form of the fh distribution determined 
from our sample of early-type galaxies is able to repro- 
duce satisfactorily the form observed in the Milky Way 
although the mean value is differs slightly.^" The dash- 
dotted curve shows the same curve after setting ^ = 0.37, 
and shows a very good agreement. This exercise supports 
the flexibility of the proposed form for describing the GC 
size distributions in galaxies of widely different mass and 
morphological type. 

■^^ This mismatch would be reduced if we had used the distance 
inferred to Virgo from the direct geometric calibration using the 
Galactic GC system 




\ [pc] 

Fig. 11. — Distribution of logf^i for Galactic GCs satisfying 
My < -7.4 mag and 0.75 < < 10 pc (66 GCs). The solid curve 
shows a Gaussian kernel density estimate with a bandwidth h = 
0.06, while the dotted curves shows its 99% confidence bands. The 
dashed line represents the distribution given in equation 1281 while 
the dash-dot-dot curve represents the same curve with fi = 0.37. 



It is worth emphasizing that we have only modelled 
the observed distribution of fh, which is the convolution 
of the intrinsic distribution with the error distribution. 
If the measurement uncertainties are large compared to 
the width of the underlying distribution, then the ob- 
served distribution might be telling us more about the 
measurement errors than the actual distribution of rh- 
However, as noted in § 12 the typical random uncertain- 
ties for our sample are of order 0'.'003 or w 0.25 pc at the 
mean distance of Virgo, with the systematic uncertain- 
ties estimated to be of the same order. The magnitude 
of the errors will thus affect the width of the observed 
peak to some extent, although the extended, high-r/j tail 
will be largely unaffected. 

To investigate the effect of random uncertainties in 
a more quantitative way, we have taken the best- 
fit function of the form of equation for VCC1316 
(M87=NGC 4486) and determined the best match, via 
least squares, of a function with this same form but con- 
volved with a Normal distribution with a — 0.25 pc in 
order to simulate the effect of the errors. The use of a sin- 
gle galaxy for this comparison ensures that we consider 
only random uncertainties: i.e., we are probing only the 
shape of the distribution, and not the absolute value of 
its location parameter. FigurelT^shows the result of this 
exercise. This figure demontrates that the inferred intrin- 
sic distribution shares the same features as the one after 
convolution, suggesting that the observed distribution is 
indeed revealing intrinsic properties. Additionally, the 
fact that J also provides a close match to the Galactic 
data, which has a completely independent error distribu- 
tion, further shows that the form of J traces the intrinsic 
properties of the fh distribution. 

Strictly speaking, the above analysis has relied on GCs 
that are brighter than the peak of the GC luminosity 
function, so our results apply in that regime only. While 
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\ [pc] 

Fig. 12. — The solid line shows the best-fit function J' to the 
distribution of VCC1316 (M87 = NGC 4486). The dotted line 
shows a function of the same form which best matches the ob- 
served data for M87 after convolution with the average random 
uncertainty in the measurement (assumed to be a Normal dis- 
tribution with (T = 0.25 pc). The convolved function is shown as 
the dashed line. 



it is reasonable to expect the distribution to apply to 
the less luminous GCs, this claim needs to be tested. 
Indeed, if we considered all Galactic GCs in Figure ITTl 
then the agreement would be somewhat worsened due 
to the presence of a larger proportion of faint, extended 
clusters. While the form of the distribution is still well 
described by J, the shape and location parameters need 
to be changed somewhat from those presented in equa- 
tion[5S| In particular, the extended high-r/j-tail becomes 
more important when the fainter GCs are included in the 
analysis, which is reflected by the larger (3i (i.e., 0.16, 
as opposed to the value of 0.12 obtained when only the 
bright GCs are considered). 

The form of the distribution presented above should 
serve are a useful constraint for models of GC forma- 
tion and numerical simulations of their dynamical evo- 
lution. Any viable picture of star formation in clusters 
should produce an observed size distribution that is con- 
sistent with the form of J', at least for GCs brighter 
than the luminosity function turnover. The fact that its 
form seems to be uniform across the wide range in lumi- 
nosity and color spanned by AC SVC S program galaxies 
strongly suggests that some quite generic mechanism de- 
termines J, while the location of the distribution is set 
by the local properties of the GC system and/or its host 
galaxy 

7. SUMMARY & CONCLUSIONS 

We have analysed the half-light radii of GCs in early- 
type galaxies using the exceptional dataset provided by 
the ACS Virgo Cluster Survey. We have presented ev- 
idence that the average half-light radii (r/j) of GCs be- 
longing to early-type galaxies increases with galactocen- 
tric distance or, alternatively, with decreasing surface 
brightness. This dependence is significantly shallower 
than that observed in our Galaxy. For the first time, we 



report an increase in (r/j) with decreasing galaxy color 
decreasing luminosity) , albeit with modest statistical sig- 
nificance. In agreement with previous observations, we 
find strong evidence that (rh) is independent of the lu- 
minosity (mass) of GCs. We also find (r^) to be ~ 17% 
smaller for red (metal-rich) than for their blue (metal- 
poor) counterparts. This trend is better explained by a 
mechanism that is intrinsic to GCs rather than to projec- 
tion effects. As such, the observations favor the picture in 
which this difference follows from the combined effects of 
mass segregation and the metallicity dependence of stel- 
lar lifetimes under the assumption of invariant average 
half-TTiass radii (Jordan 2004). 

We have discussed the predictions of two simple pic- 
tures for the origin of the of GCs: one based in the 
idea that they are determined by the average ambient 
pressure surrounding proto-GC clouds, and other that 
rfi is related to the tidal radius set by the gravitational 
field of the host galaxy. Both pictures are roughly con- 
sistent with the observed behavior, (rh) ~ L^, although 
some elements of the theoretical framework are lacking 
in both scenarios. While our level of understanding of 
GC formation is still in its infancy, it is clear that the 
observations of the rh of GCs can provide important in- 
sights into the physical mechanisms which have shaped 
their formation and evolution. 

We have studied the form of the observed distribution 
of GC half-light radii in early- type galaxies. Once cor- 
rected for a variety of dependencies on color, radius and 
surface brightness, the form of this distribution is re- 
markably consistent across the range of color and lumi- 
nosity spanned by our sample of galaxies. We introduce 
simple analytic expressions (equations 1221 and I24|l which 
satisfactorily describe the form of the fh distribution in 
our sample galaxies and in the Milky Way. Future the- 
ories of GC formation and evolution should be able to 
produce results consistent with this form when predict- 
ing the observed distribution of r^- 

Once corrected for dependencies on color and surface 
brightness, the mean half-light radii of GCs can be used 
as a standard ruler for distance estimation. Specifially, 
we find a constant value of (fh) — 2.7 ± 0.3 pc for a 
GC with color {g — z) — 1.2, in a galaxy with color 
(g — z)gai = 1.5, and at an underlying z-band surface 
brightness of /i^ = 21 mag arcsec"^. A first attempt 
at a geometric calibration of this method based on the 
Galactic GC system gives an independent distance to the 
Virgo Cluster of 16±2.3 Mpc. While these conclusions 
rest heavily on two datapoints besides the bulk of galax- 
ies in Virgo (VCC 731 and NGC 4697), the ongoing ACS 
Fornax Cluster Survey (Jordan et al., in preparation), 
which is similar to the ACSVCS but focusing on 43 early- 
type members of the Fornax Cluster, will be invaluable 
in this regard by adding many additional datapoints at 
D - 19 Mpc (Tonry et al. 2001). 

While not as efficient some other distance indicators, 
such as SBF, the use of {fh) as a distance indicator could, 
in principle, offer some advantages that could render it 
a valuable complement to other GC-based method such 
as the GC luminosity function. For instance, half-light 
radii are expected to remain roughly constant during the 
evolution of a GC. This is certainly not the case for lu- 
minosity, which will be affected by age and metallicity, 
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quantities which are not necessarily known a priori. Ad- 
ditionally, it is not necessary to go past the turnover of 
the luminosity function if enough clusters to accurately 
determine (r/j) can be accumulated up to brighter mag- 
nitudes. Assuming one could measure (fh) up to 1/4 of 
the PSF FWHM of a diffraction-limited telescope of di- 
ameter c? at a wavelength of 1/i, this technique could be 
apphed out to distances of ~ 50 Mpc {d — 6m), ~ 250 
Mpc {d = 30m) and ~ 850 Mpc {d = 100m), illustrating 
the prospect for its implementation with future observing 
facilities. In practice, the limitation will be most likely 
set by the difficulty of accurately characterizing the PSF 
over a sufficiently large field of view. 
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APPENDIX 

KINGPHOT: AN ALGORITHM FOR THE MEASUREMENT OF GC PHOTOMETRIC AND 
STRUCTURAL PARAMETERS USING HST/ACS IMAGES 

Measurements of the photometric and structural parameters of GCs in the ACSVCS, including their half-light radii, 
were carried out with code specifically written for these purposes. This code, hereafter referred to as KINGPHOT, was 
designed to perform cfhcient and accurate r^ measurements for the thousands of GCs detected in the ACS Virgo Cluster 
Survey. The primary technical challenge in the design of the code is to recover r^ for the case of marginally resolved 
objects: i.e., objects with intrinsic sizes comparable to, or smaller than, the instrumental PSF. In this appendix, we 
outline the methodology used by the code and present some simulations that illustrate its performance. 

What is observed at the detector is an integral over each pixel of the PSF-convolved object; note that we include in 
the term PSF the effects of the telescope optics and detectors. The exact nature of the physical effects that give rise 
to the observed light distribution are of no concern to us here, as long as we have an empirical way of determining 
the final distribution on the detector of the light from a celestial point source. With this definition of the PSF, V, the 
observed flux O at a given pixel in the detector is given by 

= / / dxidx2[r ®T]{xi,x2) (Al) 

J J pixel 

where T is the flux distribution of the object and (E) denotes a convolution. 

If J- can be described parametrically by a set of parameters 0, and if a model for V is available, then the parameters 
in can be simply estimated by minimization of the observed two-dimensional light distribution minus a predicted 
model. Because the family of steady-state dynamical models presented by King (1966) are known to provide an 
excellent match to the observed light profiles of most Galactic GCs, we have chosen to parameterize T by these models. 
King models are a one-parameter family indexed by the central potential, Wq, or equivalently, by the concentration, 
c = logiQ{rt/rc), where rt is the tidal radius (the point at which the density is zero) and rc is the core radius (the point 
at which the density falls to roughly half of its central value). In addition to c, there are two arbitrary scaling factors: 
the total magnitude, m, and one related to spatial scale. The latter is typically taken to the core radius, but we choose 
to work with the half-light radius r^. For King models, the ratio r^/rc is uniquely determined by the concentration, 
so the choice of either is equivalent. 

There is no analytic expression for the spatial density corresponding to the lowered Maxwellian phase space distri- 
bution assumed by King (1966), but the density is easily obtained by an integration of Poisson's equation. It would 
be prohibitively slow to compute this density profile by solving Poisson's equation during the minimization, as this 
would entail the computation of hundreds of models for each of the ^ 10^ GC candidates in the ACSVCS. Instead, a 
library of King-model density profiles was constructed, computed in a grid of 300 equally spaced points in Wq ranging 
from Wq = 1 to Wq = 12, which corresponds to the range 0.296 < c < 2.739. This library is read into memory at the 
start of KINGPHOT and a model of a given c is obtained by spline interpolation between the two closest models in 
the library. 

To represent V, we adopt the method implemented in the DAOPHOT 11 program (Stetson 1987; 1993). We refer 
the reader to Stetson (1987) for the details of the PSF representation adopted. The PSFs for the current work were 
determined empirically using moderately crowded fields in the Galactic GC 47 Tucanae as described in Jordan et al. 
(2004a). 

Given the adopted King model, the corresponding model parameters, and the PSF model, the predicted light 
distribution at the detector can be obtained by calculating the convolution in equation lAll The convolution is 
performed via a fast Fourier transform using routines from the Fastest Fourier Transform in the West library (FFTW; 
Frigo & Johnson 2005) 

Let us denote the observed data at a pixel (i, j) by dij , and the prediction of a PSF-convolved King model centered at 
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position {xo,yo) in the detector by Oy (c, m, r/i, xq, yo)- The parameters (c, to, r/i, xq, yo) are estimated by minimizing 
the value of given by 



where cr^ is the uncertainty of dij and Wij are radial weights. The value of cry used by KINGPHOT is the expected 
rms in the data image given by equation (2) of Jordan et al. (2004a). In the above summation, {|(«, j)| < rat} denotes 
all pixels within a distance rtn of (xq, j/o)- The weights Wij are set to the value of the function 5/(5 + 1/(^1^/7"^^- — 1)), 
where is the distance of the pixel's center to {xo,yo)- These radial weights are included to assist convergence as 
discussed in Stetson (1987). Initially, the parameters for all GCs were determined using rat = 4 pixels. Our simulations 
show that for extended clusters, with r/j > rj.fit/2, the measured r/j starts to suffer from biases; for larger GCs, we 
increased the value of rj-fit and redetermined the best-fit parameters. This was repeated until the condition r^ < fifit/2 
is met or rat = 15 pixel. This condition is satisfied for all GCs considered in this analysis. 

Minimizations were carried with CERN's Minuit package. The minimization was done in two steps. First, a Simplex 
(Nelder & Mead 1965) minimization algorithm was used, followed by a variable metric method with inexact line search 
(MIGRAD). In the Simplex minimization the concentration is held fixed at c = 1.5, a typical value for Galactic GCs 
(Harris 1996). The starting values for the parameters (m, xq, j/q) are obtained from SExtractor (Berlin & Arnouts 1996) 
while r/j is set to 0'.'0385, roughly the angle subtended by 3 pc at the Virgo distance. After the Simplex minimization, 
the center (xq, ya) is fixed while c is allowed to vary. Before starting the MIGRAD minimization, a grid search is done 
in (c, rh) and the minimization is started at the parameters that give the lowest this grid search. This procedure 

is done to avoid biasing the final result by the choice of fixing c in the Simplex minimization and to avoid converging 
to a local, rather than a global, minimum. After MIGRAD is run, the best fit parameters are recorded. To estimate 
the uncertainties a bootstrap procedure is used. From the oberved data dij within the fitting radius, a random set of 
dij is generated and the minimization is re-run, 50 times, with this new set. The output parameters recorded, and the 
rms of the recovered parameters is taken to represent the uncertainties. 

To test the performance of this code, it has been run on a set of simulated GCs to which an appropriate amount 
of instrumental and Poisson noise has been added. For all combinations of (c, rh) in the grid c — 0.5, 0.7, . . . , 2.1, 2.3 
and rh — 0.5, 1, 2, 4, 8, 16 pc (at the distance of Virgo), 15 objects were simulated and the best-ft parameters recorded. 
This was done for instrumental magnitudes of to = 20,21,22,23,24,25 (for a Is exposure ), wit h a fitting radius of 
rat = 4 WFC pixels. The relevant results of these simulations are summarized in Figure IA13I where we show the 
average Arh/rh (for all simulated c values). Here, Arh is the difference between the recovered and input rh, Two sets 
of error bars are shown for each point, slightly displaced. The one to the left shows the observed rms of Ar^/rh while 
the one to the right indicates the average uncertainty of Arh/rh returned by KINGPHOT. Instrumental ma gnitu des 
can be transformed into z-band ones by m^ ~ m — 0.14. Thus, only the four uppermost panels of Figure IA13I are 
relevant to the present work (due to the construction of the sample; see Equation 

We conclude on the basis of these simulations that the performance of KINGPHOT is, for observational material 
similar to the ACSVCS, largely unbiased for rh ^ 1 pc. At the same time, it is apparent that, for r?i ^ 1 pc, there 
are always systematic biases which are of order <^ 20%. It is also apparent is that a bias arises when the intrinsic rh 
approaches rat, as is revealed by the rightmost point in each panel. These bias, however, is avoided if r^, <^ 0.5rfit, so it 
is necessary that this condition is satisfied post-facto by iteratively adjusting rat as we have done (see above). Finally, 
the uncertainties returned by KINGPHOT agree very well with the observed dispersion on the Arh/rh values, and thus 
we conclude that our bootstrap procedure estimates the uncertainties properly. Needless to say, such simulations can 
provide information only on the intrinsic performance of the code, and neglect sources of error arising from inaccuracies 
in the assumed forms for J-' or V. Uncertainties arising from systematic errors of this sort can be estimated from the 
independent measurements in the two bandpasses (see §2). 





(A2) 
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Fig. A13. — Ar^j/r^ versus r/^ for a set of KINGPHOT runs on simulated GCs. The instrumental magnitude m (for a Is exposure) of 
the simulated GCs is indicated in each panel. For each simulated rfi, we plot the fractional mean difference between the measured and 
input r^, Arjj/r^j. We also show two sets of error bars for each point: the left one is the observed rms of Ar^/r;, while the right one is the 
average uncertainty computed by the code. The fitting radius was set to rflt = 4 pixels. 



